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The trace anomaly of gluodynamics encodes the breakdown of classical scale invariance due to interactions 
around the deconfinement phase transition. While it is expected that at high temperatures perturbation theory 
becomes applicable we show that current lattice calculations are far from the perturbative regime and are dom- 
inated instead by inverse even power corrections in the temperature, while the total perturbative contribution 
is estimated to be extremely small and compatible with zero within error bars. We provide an interpretation in 
terms of dimension-two gluon condensate of the dimensionally reduced theory which value agrees with a simi- 
lar analysis of power corrections from available lattice data for the renormalized Polyakov loop and the heavy 
quark-antiquark free energy in the deconfined phase of QCD [1|2) . 



Introduction. The Lagrangian of gluodynam- 
ics is conformal invariant, reflecting the absence 
of an explicit scale. The divergence of the di- 
latation current equals the trace of the improved 
energy-momentum tensor 0^ [3 and vanishes 
classically. Quantum-mechanically yields instead 
the so-called "trace anomaly" [4]. It reflects the 
breaking of scale invariance which introduces a 
single mass scale, Aqcd- The dimensionless "in- 
teraction measure" A = Tdrip/T^) = (e— 3p)/T'* 
quantifies the departure from the conformal limit 
e = 3p, which corresponds to a gas of free mass- 
less particles. At finite temperature, the energy 
density e and the pressure p enter as |5l6l7l8j , 

T^A^e-?,p^^{{Gl,f)^{Q^), (1) 

where G^i. = d^j.A^ - d^A^ - is the 

field strength tensor and I3{g) = fidg/dp = 
—bog^ + 0{g^) is the beta function, with bo — 
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lliVc/(487r^). A good knowledge of A is crucial 
to understand the deconfinement process, where 
the non perturbative (NP) nature of low energy 
QCD seems to play a prominent role [S]- In this 
contribution we analyze the highly NP behaviour 
of the trace anomaly just above the phase tran- 
sition and describe it in a way that is consistent 
with other thermal observables (see pi)! for fur- 
ther details). 

Low and high temperatures. At low tem- 
peratures A is dominated by the lightest con- 
fined states in the spectrum. In gluodynamics 
the lightest glueball mass tog ~ 1.3 GeV is much 
heavier than Tc « 270 MeV, and the pressure is 
p ^ 6"™°/^, so A ~ e"™''/^, T <^Tc, indi- 
cating an exponentially small violation of scale in- 
variance At very high temperatures one also ex- 
pects scale invariance to be restored while asymp- 
totic freedom guarantees the applicability of per- 
turbative QCD (pQCD). Actually, from the pres- 
sure to two loop one has [H] 

A= ^'^^^^~^\ o9\t^) + O{g'), T»T, (2) 

where l/g^(/i) — bolog{fi'^ / A^f^j^). It should be 
noted the ambiguity in this result, since generally 

^In full QCD for massless quarks one has a gas of weakly 
interacting massless pions and A ~ ^^/^qcd dictated 
by chiral symmetry 1111 . 
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one has both the temperature T and the (MS)- 
renormahzation scale, /i, for which one takes the 
reasonable but arbitrary choice /i ~ 2TrT — AnT. 
Higher order corrections including up to log g 
can be traced from [13^. The infrared problems 
of the perturbative expansion yield poor con- 
vergence at the lattice QCD available tempera- 
tures T < 6Tc. Perturbation theory contains 
only logarithms in the temperature, suggesting 
a mild temperature dependence. This feature is 
shared by hard thermal loops (HTL) and other 
resummation techniques of infrared divergencies 
(see e.g. |14|15] ). Actually, the value they find 
Ahtl = ± 0.5 for T > Tc is compatible with 
zero within uncertainties. Furthermore, it is not 
clear at what temperatures is the pQCD result 
dominating A. 

Thermal power corrections in gluody- 
namics. The interaction measure on the lattice 
for gluodynamics [T^ is shown in Fig. [1] and, as 
expected, is very small below Tc. It increases sud- 
denly near and above Tc by latent heat of decon- 
finement, and raises a maximum at T w 1.1 Tc. 
Then it has a gradual decrease reaching small 
values at T — 5Tc. The high value of A for 
Tc < T < (2.5 — 3)Tc corresponds to a strongly 
interacting Quark-Gluon Plasma picture. 

In previous works |l|2j we have detected the 
presence of inverse power corrections in other 
thermal observables. In Fig.[l]we plot (e — 3p)/r'* 
as a function of l/T^ (in units of Tc) exhibiting 
an obvious straight line behaviour in the region 
slightly above the critical temperature. 
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and corresponding to a "power correction" in 
temperature. A fit of the lattice data {N^ x Nt^ = 
32^ X 8) for 1.13rc < T < A.54Tc yields atra = 
-0.02(4), 6tra = 3.46(13), xVDOF = 0.35. 
Power corrections also appear in e and p, just by 
applying the standard thermodynamic relations. 
Note that if the power correction was absent the 
transition towards the free gas value would take 
over much faster indicating a weakly interacting 
quark-gluon plasma in the neighbourhood of Tc. 

Pisarski [IT] has suggested to interpolate be- 
tween T = Tc where p ^ for non-interacting 



and heavy glueballs and the free massless gluon 
gas value p ^ T'^ at T :s> Tc. If we choose 
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T>Tc, (4) 



with n arbitrary we get 
n{N^ - l)7r2 ^T 



A = 



45 



T 



T>Tc 



(5) 



which for n ~ 2 corresponds to take = and 
b — 3.51 in Eq. ([3]), in excellent agreement with 
the fit to the lattice data. This thermodynamic 
consistency does not explain however why there 
is a power correction with n = 2. A fit of Eq. ([5]) 
to the data of Fig. [1] for the same range as in 
Eq. ©, yields n = 1.97(5) with xV^OF = 0.62, 
indicating the robustness of the power. 

The lattice behaviour of A clearly contra- 
dicts perturbation theory [12ll3j and resumma- 
tions thereof |14I15| explaining why they have fla- 
grantly failed to describe the data of the free en- 
ergy below 3 Tc. Our discussion above shows that 
these approaches would yield a powerless contri- 
bution. Apt, which should ultimately be identi- 
fied with the almost constant and vanishing atra 
of Eq. ([3]) rather than with the full result from the 
lattice |16| . Actually, the maximum lattice tem- 
perature T — 5Tc may still be far from the pQCD 
estimate since the power correction provides the 
bulk of the full result at this temperature. For 
TVc = 3 the 0{g^) correction to Eq. © corre- 
sponds to multiply it by (1 — 6g{^)/'K) [13| which 
becomes small for g{fi) ^ n/6 or » IO^^Aqcd- 
This delayed onset of pQCD is not new and hap- 
pens in the study of exclusive processes at high 
energies (see e.g. Ref. [H]). 

Dimension two gluon condensate. Power 
corrections are unmistakable high energy traces 
of NP low energy effects. Within QCD sum 
rules power corrections at high are usually 
related to local condensates as suggested by the 
OPE. The gluon condensate (G^) = g^(iGl^f) 
describes the anomalous (and not spontaneous) 
breaking of scale invariance, and hence is not an 
order parameter of the phase transition. Actually, 
the order parameter is the vacuum expectation 
value of the Polyakov loop L{T) which signals the 
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Figure 1. Trace anomaly (e — 3p)/T^ as a function of T (left) and l/T^ (right) (in units of Tc). Lattice data are 
from [in] for x Nr = 16^ x 4 and 32^ x 8. The fits use Eq. ^ with otra and 6tra adjustable constants. 



breaking of the 'Z{Nc) discrete symmetry of gluo- 
dynamics as well as the deconfinement transition. 
A dimension two gluon condensate naturally ap- 
pears from a computation of L{T) which in the 
static gauge, doAo{x, sq) = 0, in a Gaussian- like, 
large Nc motivated, approximation gives 



1 

Nr. 



— tree^f^«(^)/^ 



exp 



ANrT^ 



(6) 



valid up to 0{g^) in pQCD. is the gluon field 
in the (Euclidean) time direction. 

The dimension two gluon condensate g^(^oa) 
is obtained from the gluon propagator of the di- 
mensionally reduced theory, -Dooj by taking the 
coincidence limit. The perturbative propagator 
DP (k) = l/(k2 + ml) + 0{g^), being mo ^ T 
the Debye mass, leads to the known perturbative 
result [19j and fails to reproduce lattice data be- 
low 6 [ini- A NP model is proposed in Ref. [T] 
to describe the lattice data of the Polyakov loop, 
and it consists in a new piece in the gluon prop- 
agator driven by a positive mass dimension pa- 
rameter: 



Doo(k) = DP(k)+Do^P(k) 



(7) 



where ^'□□^(k) — m^/iyi^ + ■ This ansatz 
parallels a zero temperature one [21j . where the 



dimension two condensate provides the short- 
distance NP physics of QCD and at zero temper- 
ature this contribution yields the well known NP 
linear term in the qq potential. A justification of 
Eq. ([7]) based on Schwinger-Dyson methods has 
been given [22j . The Gaussian approximation has 
also been used in Ref. [2] to compute the singlet 
free energy of a heavy qq pair [20|23j , through the 
correlation function of Polyakov loops. 

Non perturbative contribution to the 
Trace Anomaly. The model of Eq. ^ can eas- 
ily be used to compute the trace anomaly Eq. ([1]) 
in gluodynamics. Assuming the leading NP con- 
tribution to be encoded in the ^o,a field and tak- 
ing A,,a = yields [24] 



2\NP 



(8) 



Note that the NP model is formulated in the di- 
mensionally reduced theory, so the gluon fields 
are static. This formula produces the thermal 
power behaviour of Eq. ([3]) with 



6traT2 = -(3/3(5)/.g)m2,(A2,J 



NP 



(9) 



where mo = tiid/T. If we consider the perturba- 
tive value of the beta function (3{g) ^ g'^ + 0{g^), 
the r.h.s. of Eq. ([9]) shows a factor g^ in ad- 
dition to the dimension two gluon condensate 
5^(^0a)^^- ^'^ the fit of the trace anomaly 
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data is sensitive to the value of the smooth T- 
dependent 5, without jeopardizing the power cor- 
rection. When we consider the perturbative value 
gp up to 2-loops, we get from the fit of the trace 
anomaly a value of g'^{A^ a)^^ which is a factor 
1.5 smaller than from other observables. This dis- 
agreement could be partly explained on the basis 
of certain ambiguity of g in the NP regime. A bet- 
ter fit of the Polyakov loop and heavy quark free 
energy lattice data in the regime Tc < T < ATc 
is obtained for a slightly smaller g than gp, i.e. 
g = 1.26 — 1.46 2J. Taking this value we get 
from Eq. ^ 5^(^o,a)^^ = (2-86 ± 0.24 Tc)^, 
a better overall agreement, see Table [T] From 
Renormalization Group requirements it is possi- 
ble to observe that the consistency of the model 
is warrantied if the coupling constant as{fJ.) = 
5^(/i)/(47r) has a behaviour 1//^^ at low enough 
temperature [lOj . This behaviour is just what is 
obtained within the Analitic Perturbation Theory 
formalism, after extracting the Landau pole [25j , 
with the corresponding decrease of the perturba- 
tive value of as (fJ,) . This decrease could approxi- 
mately explain the best-fit value of g. 

While there might exist other explanations to 
the observed thermal power corrections our re- 
sults of Table [1] suggest an unified and coherent 
description of observables in the non perturbative 
regime of the deconfined phase (sQGP) in terms 
of the dimension two gluon condensate [lOj . 



Observable I g^(Ag^)^^^ 
Polyakov loop llj (3.22 ± 0.07 T^)^ 

Heavy qq free energy [2] (3.33 ± 0.19 tJ}^ 

Trace Anomaly | (2.86 ± 0.24 tJ)^ 

Table 1 

Values (in units of Tc) of the dimension two gluon 
condensate from a fits in the deconfined phase of gluo- 
dynamics: Polyakov loop, singlet free energy of heavy 
quark-antiquark and trace anomaly. Lattice data are 
with Nt = 8. The error reflects an uncertainty in 
the coupling constant g = 1.26 — 1.46, being the 
highest value the perturbative gp up to 2-loops at 
T = 2Tc. The critical temperature in gluodynamics 
is Tc = 270 ± 2MeV [20] . 
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